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In this paper we examine in detail the implementation, with its associated difficulties, of the 
Killing conditions and gauge fixing into the variational principle formulation of Bianchi-Type cos- 
mologies. We address problems raised in the literature concerning the Lagrangian and the Hamilto- 
nian formulations: We prove their equivalence, make clear the role of the Homogeneity Preserving 
Diffeomorphisms in the phase space approach, and show that the number of physical degrees of 
freedom is the same in the Hamiltonian and Lagrangian formulations. Residual gauge transforma- 
tions play an important role in our approach, and we suggest that Poincare transformations for 
special relativistic systems can be understood as residual gauge transformations. In Appendices, we 
give the general computation of the equations of motion and the Lagrangian for any Bianchi-Type 
vacuum metric and for spatially homogeneous Maxwell fields in a nondynamical background (with 
zero currents). We also illustrate our counting of degrees of freedom in an Appendix. 
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I. INTRODUCTION 

A spatially homogeneous cosmological model is a mani- 
fold M with a Lorentzian metric tensor g invariant under 
a group of isometries whose three-dimensional, spacelike, 
invariant hypersurfaces foliate A4. In the models we will 
treat, this group is generated by three spacelike vector 
fields Ka, which span a Lie algebra defined by their com- 
mutation relations: 

[K„,Kf,] ==C^,K, . (1.1) 

The invariance of the metric is expressed by the vanish- 
ing of its Lie derivatives with respect to these vectors: 
'CkbS = 0. The structure constants are antisymmetric in 
their lower indices and obey the Jacobi relation: 

'^ab — ~'^ba i '^be'^cd + '^ce^db + ^de^bc " ■ U-^J 

The Bianchi classification of these algebras into nine 
Types (see |l]) is the source of the term Bianchi-Type 
cosmology. In three dimensions the Jacobi relation is 
equivalent to 

C^aQc = . (1.3) 



The algebras with vanishing C"^ are called Class A and 
those with nonvanishing C^^ are called Class B . 

In a suitable basis, g can be represented by components 
gfj_i, which depend only on a single variable, cosmic time t. 
The Einstein field equations are (coupled, non-linear) or- 
dinary equations. (The choice of this basis is part of this 
paper.) The Einstein-Hilbert Lagrangian density for the 
metric may be computed in terms of these components 
and their ^-derivatives. It is well known that in Class A 
models this reduced Lagrangian does correctly reproduce 
the field equations, but it does not necessarily do so in 
Class B models because of a spatial divergence that 
is automatically zero only in Class A 

Bianchi-Type cosmology is still a subject of debate 
and — we think — some misunderstanding. It has been 
claimed, surprisingly, that the number of degrees of free- 
dom in the Lagrangian and Hamiltonian formulations of 
these models do not agree Q. This result would mean 
that the two formulations are not physically equivalent. 
Another issue, related to this one, poses the question 
as to what is to be considered a gauge transformation, 
and therefore a redundancy in the physical description in 
these models. 
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In this paper, we will address three subjects: the La- 
grangian and Hamiltonian formulations of Bianchi-Type 
models (dimensional reduction), the concept of gauge 
freedom for these models, and the fact that the number 
of degrees of freedom in the Lagrangian and Hamiltonian 
formulations do agree. We do not solve the first subject, 
which remains a problem for Class B models, although 
we clarify some points that are also relevant for Class 
A models. From our analysis we give an answer to the 
second and third problems. We show that the gauge free- 
dom is dictated by the diffeomorphism invariance of the 
original theory and leads to the physical equivalence of 
the Lagrangian and Hamiltonian formulations. 

We treat homogeneous models in general relativity as 
constrained dynamical systems. In particular we clar- 
ify the two aspects that merge when we move from the 
superspace of all metrics to the minisuperspace of ho- 
mogeneous metrics. On the one hand there is the di- 
mensional reduction from 3 + 1 spacetime dimensions to 
one time dimension. This dimensional reduction comes 
from the symmetry associated with the Lie algebra of 
the Killing vectors. On the other hand there is the elim- 
ination of gauge degrees of freedom. These two types of 
reduction present different types of problems, which we 
analyze, but both are necessary in order to reduce the 
original gauge group of four-diffeomorphism invariance to 
the gauge group of time reparameterization invariance. 

This reduction procedure will be undertaken in four 
steps. As an outcome of our analysis, a) we will clarify 
some points concerning the Lagrangian and Hamiltonian 
reduced formulations; b) we will be able to point out 
clearly what is the gauge group for these models; and 
c) we will show that the correct reduction procedures in 
both the Lagrangian and Hamiltonian formulations lead 
to a number of degrees of freedom which is always the 
same in either formulation. The first step, in Section 
consists of adapting time coordinates to the symme- 
try group. The second step, in Section is to adapt 
spatial coordinates to the group by adopting a gauge in 
which the shift vector depends only on time. The third 
step, in Section IV, is to adopt a gauge in which the shift 
vector vanishes. Finally, the fourth step, in Section is 
to eliminate residual three-diffeomorphisms. After each 
step we examine the status of what remains of the gauge 
group. In Section VI we examine the loss of constraints 



ground metric. Appendix |c| reviews the general litera- 
ture on the relationship between setting gauge conditions 
and variational principles. It gives justification to some 
results used in Section IV and Section VII. Appendix 



caused by the gauge fixing. We perform, in Section VH 



the Hamiltonian analysis, which is presented in a sum- 
mary form because it is parallel to the Lagrangian anal- 
ysis. Section VIH is devoted to conclusions. 

In Appendix ^ we present explicit equations for spa- 
tially homogeneous metrics in full detail, including gen- 
eral lapse function and shift vector, for any Bianchi Type, 
where coordinates have been chosen so that all variables 
depend only on time. The purpose is to provide concrete 
examples for our formalism. A simpler example, which 
also can serve to illustrate our methods, is provided in 
Appendix ^ It is the case of a spatially homogeneous 
electromagnetic field in a spatially homogeneous back- 



Pl illustrates our counting of degrees of freedom in the 
Bianchi Type I case and shows how spacetime rigid sym- 
metries may be considered as residual gauge transforma- 
tions. 



II. TIME 

We start with the general setting for Bianchi models: 
Call T the set, each element of which is a metric ten- 
sor g and three vector fields Kq whose orbits generate 
three-dimensional hypersurfaces which foliate the four- 
dimensional manifold A4; the Lie derivatives of g with 
respect to the Kq vanish. The invariance of g and even 
the definition of the Ka need only be locally defined for 
much of what we do; if the vector fields are globally de- 
fined, we speak of global homogeneous cosmologies [|| . 

The diffeomorphisms on M can be viewed in an active 
sense, mapping points onto other points. In that case, 
an element (metric plus three vector fields) of T will in 
general be mapped to another element in T. If there is 
a covering of Ai by coordinate patches, each diffeomor- 
phism may also be viewed in a passive sense as a col- 
lection of coordinate transformations. It is this passive 
sense which is closer to the physical principle which is 
a motivation for General Relativity, that physics should 
not depend on coordinates. It is, in fact, somewhat easier 
to adopt the language of a metric being determined by its 
components in a coordinate patch; a diffeomorphism is 
then a transformation which preserves the metric g and 
the vectors Ka but in general changes their components 
gfiu^Kct. Nevertheless, we shall adopt the active view 
when convenient. (Greek indices range over 0,1,2,3, with 
the coordinate x'^ being the time t. Latin indices will 
range over 1,2,3.) 

Given a coordinate neighborhood A/", we consider all 
metric tensors, defined by their components g^j/ in A/", 
which are invariant under isometrics defined by the 
and in which the are spacelike. The metric compo- 
nents satisfy the Killing equation: 

(2.1) 

where comma denotes partial differentiation. We can 
consider T as being the collection of all such tensor com- 
ponents, each element of T being the collection of metric 
and Killing vector components, the coordinate system 
being understood. Any of the elements of T foliates Ai 
by three-dimensional space-like homogeneous hypersur- 
faces, namely the integral surfaces of the Killing vector 
fields. 

A change of coordinates — a diffeomorphism — will 
in general change the form of the metric components. In 
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that sense, the diffeomorphism group Diff (Al) on A4 is 
reahzed as a group D acting on T. The physics doesn't 
change, of course, under a coordinate transformation, 
and symboHcally we can write 

Physics {of a Bianchi Type) — T/D . 

These conditions need only hold locally for much of 
what we will be doing. In that sense, we really are work- 
ing with a Lie algebra of infinitesimal isometries rather 
than a Lie group, though we will continue to speak of 
the isometries as a symmetry group. When we need to 
require global homogeneity, we will clearly specify so. 

The existence of these isometries of the metric can be 
used to simplify greatly the equations of motion, which 
we take to be the vacuum Einstein equations. At this 
point, the full gauge group D of General Relativity, gen- 
erated by DifFAl, is still operating in T. We will be look- 
ing for a reduced Lagrangian, with the Killing conditions 
built in, capable of describing these Bianchi models. We 
do so in four steps: adapting time, adapting space, fix- 
ing the gauge by requiring the shift vector (to be defined 
later) to vanish, and fixing the residual gauge freedom. 

We now proceed to the first step, which will partially 
fix the gauge (that is fixing coordinates) by concentrat- 
ing on the three-dimensional homogeneous hypersurfaces. 
Suppose that the homogeneous hypersurfaces happen to 
coincide with the hypersurfaces t = const in the coor- 
dinates defining an element of T, {g^,y,Ka)\ then we 
introduce an equivalence relation by saying that another 
element of T is equivalent to this one if it is produced by 
an element of D, that is by a diffeomorphism. The set 
of these equivalence classes is Tt, and the full group D 
has been reduced to a smaller group Dt, namely those 
diffeomorphisms which preserve the condition that the 
homogeneous hypersurfaces are defined hy t — const. 

In this first step, we have made a choice of spacetime 
local coordinates such that the surfaces Sf of constant 
time coincide with the foliation defined by the Killing 
vectors. The time coordinate is a function whose curl 
vanishes under Lie differentiation by any of the Killing 
vectors: 



Cm At = 



(2.2) 



Our Killing vectors now take the general form (where 
9, = d/dx') 



K„ 



K^{t,x)d^ 



(2.3) 



Notice the possible time dependence of Kq. (We denote 
the four spacetime coordinates t, x^ by t, x when they are 
used as the arguments of a function.) 



A. Gauge group after the first step 

In these spacetime coordinates {t, a;*}, a general in- 
finitesimal reparameterization (coordinate transforma- 
tion) is generated by a vector field e'^(x)c?^: 



t + St 



t + e" 

x' + 



(2.4) 



To preserve equation ( p. 2D , that is, to generate a foliation- 
preserving diffeomorphism, requires 

e° = . (2.5) 



B. Expressing the homogeneous metric 

Recall t he c ommutation relations of the Killing vectors 
equation (1.1). At this point it is customary to introduce 
a set of three independent, right-invariant (under the Lie 
algebra) vector fields Ya = Yait, x)di, which satisfy 



[Ka,Yb] = . 



(2.6) 



We take them to be tangent to the homogeneous hyper- 
surfaces. They also define a Lie algebra and can be taken 
such that 



(2.7) 



These commutation relations are ensured if we take Y^ 
to coincide with at an arbitrary point in Et and define 
them at other points in Et by equation (^) |^ . 

In each we can define a basis of 1-forms, 
a;" — LLjf{t, x)dx\ dual to these left-invariant vectors: 
uj"-(Yb) — 6^. The Lie algebra property becomes (da is 
differentiation with respect to the space variables) 



(2.8) 



The metric can be written, using the anholonomic ba- 
sis {dt, uj'^}, as 



g = -N^dt^ + gab (N^'dt + w") {N^'dt + J" 



(2.9) 



where iV is the lapse function and the shift vector 
variables in this basis. The Killing conditions on g give 
( ■ denotes djdt) 



Ka(5fcc) = ■^=4> Qab = gah{t) 



KaiN) N ^ Nit) 



CKAN''dt + uj'')=0 

^ = [N''Yb, Ka] 



^Ka{N'')Yb 



(2.10a) 



(2.10b) 



(2.10c) 



The first two results are clear: no spatial dependence for 
N and gab- The third result links the time dependence 
of the Killing vectors to the spatial dependence of the 
shift variables N"'. It relates dimensional reduction and 
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gauge-fixing reduction. Since the shift variables are com- 
pletely arbitrary because of the gauge freedom, n ame ly 
diffeomorphisms satisfying e° = (see equation (2.5)), 
we may consider this third relation as giving the time 
evolution of the Killing vector components for a given 
set of the shift variables. This shows in particular that if 
for an "initial" time, say t — 0, the Killing conditions are 
satisfied, they will be satisfied in future times, the third 
relation providing us with the form which the Killing vec- 
tors take. 

At this point we could proceed to reduce the La- 
grangian (which we will describe below). But reduction 
within the Lagrangian and reduction in the equations of 
motion are procedures that in general do not commute, 
as we could immediately verify. In the next section a 
thorough study of this non-commutativity will be given. 

The fact that the shift variables depend arbitrarily on 
the space coordinates shows that we must further reduce 
the gauge freedom in order to end up with the reduced 
gauge group of time-reparameterization invariance. To 
do so, the simplest way is to require the shift variables 
to be spatially constant. According to the third relation 
above, this is equivalent to requiring time independence 
of the Killing vectors. This is our second step in the 
gauge fixing procedure, in the next section, and is a case 
of adapting spatial coordinates to the Killing structure. 



III. SPACE 



We introduce the partial gauge-fixing conditions 



N^, = 



(3.1) 



We know from the previous analysis, equation (2.10c) 
that the Killing vectors become time independent, 

Ka = K[x)d, , 



and we can also choose the right invariant vectors Yq 
and the 1-forms a;° to be so. It was realized in [|lO| that 
the shift vector N = N°-{t)Ya is a one-parameter family 
of inner automorphisms of the right-invariant Lie alge- 
bra (generators Yf,). More details are given in Section 
Now the metric is written as 



-gabit) iN-it)dt + uj''){N\t)dt- 



(3.2) 



with 



dtjj 



1 



and the Killing conditions are built in. 

We proceed to study the non-commutativity of the op- 
erations of introducing the conditions gat.i — A^,i = = 
into the equations of motion (that is, to look for a 
restricted set of solutions of the Einstein equations) or 



directly into the Lagrangian. The calculation given in 
Appendix ^ shows explicit examples of the ideas we dis- 
cuss in this Section. In the anholonomic basis {dt, w"}, 
the Einstein-Hilbert Lagrangian density, once a total di- 
vergence is discarded, can be written as (our notations 
follow 0) 



C = \u;\^N{^R + KabK^' - K') = \lo\C 



(3.3) 



where = det(u;f), g — del{gab), N ^ Q is the lapse 
function, '^R is the three-metric curvature, Kab the ex- 
trinsic curvature 



1 



Kab = ^(5afc - Na\b 



N, 



b\a) 



(3.4) 



and K = e°'^Kab its trace (| denotes the spatial covariant 
derivative and e"'' is the matrix inverse of gab)- Notice 
that C is second order in the three-space derivatives. 

One could also add other fields in the Lagrangian 
and implement the Killing conditions on them. For 
instance, in the Einstein-Maxwell case, the one-form 
gauge field A = A^dt + Aa^jj"" must satisfy Ao.i — 
Aa,i = (see Appendix ^). If we label as X the 
generic variable gab, N, N'^ , Aq, Aa, then £ is a fimc- 
tion £{X, X ,YbX,YaYbX), where the last argument in 
£ has a > b. Notice that YaYbX contains not only sec- 
ond order but also first order spatial derivatives of X . 

Let us write the Euler-Lagrange equations using the 
variables just displayed. The variation 5C is 



6C 



dC 



dX 



6X 



dC 

" dX 
dC 



6X 



dC 



i>b 



dYaYbX 



dYaX 
5{YaYbX) 



SiYaX) 



Integration by parts, dropping the boundary terms, and 
repeated use of the relation 



Y,(|c.|) + |c.|a,y,^ = |c.|c,^, 



(3.5) 



yields the following form for the Euler-Lagrange equa- 
tions: 

5C , ,f dC ^ dC ^„ dC , 



-Y^{Yb + Ct){Ya + C:,){ 



dC 



i>b 



dYaYbX 



(3.6) 



Due to the particular structure of equation (3.3), the 
third and fourth order spacetime derivatives in equation 
( |3.6D will cancel, but this fact does not alter our discus- 
sion. 

If we define the reduced Lagrangian 
Cr{X, X) = C{X, X, YaX = 0, YaYbX = 0) , 
we end up with 
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\5X 



dC . 
dC . 



a>b 



(3.7) 



equation (3^) displays the non-commutativity between 
the two procedures: Reduction of the Lagrangian 
through the Killing conditions or reduction of the equa- 
tions of motion. As we have mentioned before, this prob- 
lem was identified long ago. The distinction between 
Class A and Class B Bianchi models is crucial in this 
respect, for Class A is characterized by the vanishing of 
the trace C^f, [|. 

Summing up, we see that the implementation of the 
Killing conditions within the Lagrangian produces no 
harm, that is, commutativity holds in equation (3.7), for 
Class A models. For Class B models the equations of 
motion derived from Cn are wrong when there is a con- 
tribution different from zero coming from the last two 
pieces in equation ( |3.7| ). One may wonder whether this 
means that the situation is hopeless if we want to have a 
Lagrangian formulation for Class B models. Perhaps in 
general, but in some special cases there may be ad hoc 
solutions. In Appendix ^ we show an example of such a 
case within homogeneous Maxwell theory. 



A. Gauge transformations after the second step 

At this stage we have time-independent Killing vec- 
tors. Let E = e^c?^ = e°c?o + e'^Y^di be a generator of a 
diffeomorphi sm w hich preserves this requirement as well 
as equation (^.S]): 



[E,KJ= [E,KJ =0 , e" = eO(0 , 



(3.8) 



where ' = d/dt. The general form for the spatial compo- 
nents of E is given by 



(3.9) 



with 1^" (t) and (x) arbitrary functions of time and spa- 
tial coordinates respectively. 

Suppose El and Ei are two generators satisfying equa- 
tion 



[Ei,K,] = [E2,K, 







By the Jacobi identity it is clear that [[Ei, E2], Ka] ~ 0, 
but it is not necessarily the case that [[Ei, E2] ', Ka] — 0. 
A closure process which is more general than commuta- 
tion applies for these generators. The infinitesimal dif- 
feomorphism generated by Ei changes the Killing vec- 
tor fields (in the active view of diffcomorphisms), but it 
also changes the invariant basis vector fields and the sec- 
ond generator E2 as well. These changes, particularly 



changes in the in varia nt basis vectors Yq , mean that the 
form of equation ( |3.9| ) cannot be expected to be invariant 
during the process of commutation. We will not pursue 
this matter further, for in the next step, where the shift 
vector is set to zero, the form of the gauge group once 
again becomes straightforward. 

Notice that for Cla ss A models, the first term on 
the right in equation (3.9) is a N oeth er symmetry for 
Cr. Our Lagrangian C, equation ( |3.3| ), differs from the 
Einstein-Hilbert Lagrangian Ceh by a divergence 



C = C 



EH 



\a ■ 



General reparameterization invariance under a diffeomor- 
phism generated by v e^^{x)d^ = ip''{t)Ya produces 
the functional variation of Ceh, 

KCeh = df.ie^'CEH) ■ 

Therefore, 

6^£ = d^ie^CEH) + <5v {H^iYam'' + C^bm")) . 

(3.10) 

Let us now reduce these expressions, for Class A mod- 
els, introducing our second step gauge fixing, diX — 
dijX = 0, where X = gab, N, iV", ^0, Aa, . . . . The left 
side of equation (|3.10 ) becomes 



|w|(5v'C 



R ■ 



The first term on the right in equation (3.10) is 



(aiX=8; ■x=o) 







where we have used equation ( |3.5| ) in the l ast step. The 
second term on the right in equation (3.10) becomes 



= \.\YASAV-9-^n),,^^__,^^^__„,=o. 

We end up with SvCr = 0; this is a gauge Noether sym- 
metry provided by the reduced theory. 



The second term on the right in equation (3.9) is not 
retrievable as a gauge symmetry from the reduced La- 
grangian formulation (except when it is a copy, with con- 
stants if"" , of the first term). However, we must still pro- 
vide a gauge fixing for it. In step four of our gauge fixing 
procedure we will deal with it. It is obvious from equa- 
tion (3.9) that the process of reducing the gauge group 
is not yet finished; we now continue this task. 
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IV. ZERO SHIFT 

The remaining diffeomorphism invariance allows us to 
perform a new partial gauge- fixing by introducing a set of 
constraints that helps to eliminate some of the arbitrari- 
ness that exists in our equations of motion. To simplify 
the analysis we consider the case in which the metric is 
the only field in £. Should other gauge fields be present, 
the dynamics could have additional arbitrariness. 

Notice that the Lagrangian is independent of N and 
N"-. This fact implies that the dynamical evolution 
vector field operator in configuration-velocity space has 
| p2| , p^ a term of the type (see Appendix ^ for clarifica- 
tion of this point) 



that, besides keeping the three- foliation, preserve the 
conditions N'^ = 0. This requirement results in 



(,0(M)^-f.'^(M)^) , (4.1) 



with rp tT]"" arbitrary functions. The important conse- 
quence of this arbitrariness is that to fix the dynamics 
we must fix the values of N and N"- . Here we are only 
interested in a partial gauge-fixing: We only fix the shift 
functions. The simplest way of doing so is by introducing 
the gauge- fixing constraints N°' = 0. It is always possible 
to pass from an initial configuration with N°- to a fi- 
nal configuration with iV" = by usin g dif feomorphisms 
satisfying the restrictions of equation ( |3.9| ); the geomet- 
ric picture of this transformation is to make the curves 
in M tangent to the normal vector of the Et coincident 
with the curves generated by dt- 

This is the Lagrangian version of the partial gauge fix- 
ing. Stability of the new constraints under time evolution 
yields the new constraints N'^ = 0; then requiring stabil- 
ity again will make the arbitrary functions ry" = 0. 

To complete the picture, let us go back to the se cond 
gauge fixing step, equation (B.l). Consider equation (4.1) 
again. Stabilization of TV^ — implies iV"^ = 0. The sta- 
bilization of Nl — yields 77" — 0. By the same token, 
rfj vanishes as a consequence of the relation, which is also 
a gauge fixing, N^i — 0, obtained in the first gauge fixing 
step. Therefore, from the point of view of the reduc- 
tion of the equations of motion, the arbitrariness in the 
dynamical evolution vector field operator is described by 
four functions 77" (i), ri°-{t). The gauge freedom associated 
with this arbitrarin ess i s given by €^(t), f°'(t) in equation 
(2.5) and equation ( |3.9| ). Notice that the second term in 
equation (^|^) is unretrievable from the reduced dynam- 
ics. This is another way to verify the limitations of the 
reduced formalism under the Killing conditions already 
pointed out in the last paragraphs of Section 111 . 



A. Gauge group after the third step 

Once the shift functions have been set to zero, the 
gauge group has been greatly reduced, for the only re- 
maining diffeomorphisms still available to us are those 



= e" 



(4.2) 



The remaining diffeomorphisms therefore are such that 
e° only depends on t (see equation 2.5), whereas the e" 
depe nd only on the three-space coordinates (see equa- 



tion O). A nice picture emerges: The remaining gauge 
group has been factorized into two commuting subgroups, 
the group of time reparameterizations and the group of 
three-space diffeomorphisms. 

Properly speaking, only the first group is still a gauge 
group; it is directly associated with the freedom that is 
left in the Lagrangian evolution operator and that is dis- 
played in the arbitrariness of rf in equation ( |4.1| ) (77" are 
zero in our particular gauge- fixing) . 

The group of three-space diffeomorphisms is not a 
gauge group because it does not have room for what is 
characteristic of the gauge freedom: to connect the mem- 
bers of a family of field configurations, all of which are 
solutions of the equations of motion, that share the same 
set of initial- value conditions. This group must be under- 
stood as describing a redundancy in the space of initial- 
conditions for our theories. We call this group a residual 
gauge group. A parallel case in electromagnetism is the 
residual gauge symmetry that is left after the introduc- 
tion of the Lorentz gauge df^A^^ — 0: The transformation 
Af^ — > -I- 9^A is a residual gauge symmetry if DA = 0. 
In any case, one must take into account that the gauge 
fixing procedure is only finished when we have c omp letely 
removed these residual gauge transformations p3| . Fur- 
ther comments on this important point are made in the 
last section of this paper. 

We finish this section by observing that the comments 
raised in the previous section still apply here: Only the 
diffeomorphisms of the form e* — iJ'^Y^, with con- 
stant are obtainable as gauge Noether symmetries from 
the reduced Lagrangian in Class A models. In the lan- 
guage introduced in the next section, these vector fields 
B^Y^ define the inner automorphisms for the right in- 
variant Lie algebra (generators Ya) 



V. RESIDUAL THREE-DIFFEOMORPHISMS 

Now we are ready to perform the fourth step of the 
gauge fixing. Notice that the three-metric, gab{t)u)°'u)^ on 
the surfaces is also Killing with respect to the vectors 
Ka. We are going to fix, in principle, the residual three- 
diffeomorphism group by working from now on with this 
basis of one-forms {w"}. In the active view, a general 
three-diffeomorphism will drive Ka K'q, a;'". 
The transformed metric gab{t)u3'°'u}''^ = g' ^(^{t, x)uj°-u)'' 
will be Killing with respect to the new set K'a and will 
in general no longer be Killing with respect to the original 
vectors Ka. In other words, will in general acquire 
a dependence on the spatial coordinates. However, there 
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is another possibility to explore: the case when the new 
set K'a happens to belong to the Lie algebra generated 
by the original Killing vectors Ka- This is equivalent to 
saying that 



(5.1) 



with constant. These are the homogeneity- 

preserving difFeomorphisms (HPD), which were intro- 
duced by Ashtekar and Samuel 

Generators v for the HPD are better expressed in the 
basis of the invariant vectors Y^, v = /°(x)Ya. The 
condition for an HPD is (£ is Lie derivative here) 



with a constant matrix. This is equivalent to 



(5.2) 



(5.3) 



with a constant matrix. The left-invariant Lie alge- 
bra (generators K^) and the right-invariant Lie algebra 
(generators Yq) possess the same automorphisms. If the 
point xq in S is the point where Y^ coincides with Ka, 
then - Ai+r{xo)C',,. 
Equation (^.3[) is 



(K,r)Ye = -[v,KJ - A^K, 
^ {KaD = (K, . u:^)Ai 

This last equation is equivalent to 

df - (Kfc • u;^)A^K" , 



(5.4) 



(5.5) 



where K° are the dual forms to K^. Notice that (Kb •o;'^) 
is the adjoint representation of the Lie group expressed 
in the local patch we are working in. Because dK*^ = 
iC^fcK'^ A K*- and d{Kd ■ w^) = C^"d(Ke • u)'=)K^ the 



local integrability conditions for equation (5.5), 

d((Kb.c^^)^^K'^) =0 , (5.6) 

read 



(5.7) 



which are the conditions for Lie algebra automorphisms. 
Trivial solutions are A^; = and A^; = A^C^^, with A"" 
constant. The solutions for the first case are 



V = B^Y, , 
with B'^ constant, and for the second, 

V = A^Kc . 



(5.8) 



(5.9) 



The first (second) solutions describe the inner automor- 
phisms for the right-invariant (left-invariant) Lie algebra. 



In finite form , the in equation (5.1), (M = e , B of 



equation 5.2), satisfy 



(M-i)^(M- 



ah 



(5.10) 



The integrability conditions equation (5.7) only guar- 
antee the local existence of the HPD. But we need these 



HPD to be global diffeomorphisms if equation ( 5.10 ) is 
to hold everywhere. Both the Lie algebra structure and 
the topology of the surfaces of homogeneity play a role in 
the determination of the existing HPD. In particular, for 
globally homogeneous cosmologies with simply connected 
surfaces of homogeneity, every solution of equation (5.7) 
defines a generator of a HPD. The relevant role of spatial 
topology was first emphasized in and treated in great 
detail in JT6|-p^. The degrees of freedom for Class A 
globally homogeneous Bianchi models are also discussed 
in this reference; we will return to the question of degrees 
of freedom later. 



A. Gauge group after the fourth step 

Once we have completely eliminated the residual three- 
diffeomorphism gauge invariance, the remaining gauge 
group is that of time reparameterizations. The lapse 
variable transforms, under the time reparameterization 
gauge group, as a scalar density: 

5^a(t)N = e°N ^i^N , (5.11) 

and the three-metric variables gah as a scalar: 

^e0(t)9ah = e°3ah • (5.12) 



VI. CONSEQUENCES OF iV" = 



We have seen in Section III the non-commutativity 
between the operations of inserting the Killing condi- 
tions into the equations of motion or directly into the 
Lagrangian. Here we are going to consider another kind 
of non-commutativity, the one coming from plugging the 
gauge fixing iV° — into the Lagrangian. To make things 
clear , we only consider here the reduction of C in equation 
( |3.3|) to a reduced Lagrangian £gf, prior to the introduc- 
tion of the Killing conditions in it (GF stands for gauge 
fixing) . 

The important question is whether this new La- 
grangian £gf is going to reproduce the same equations of 
motion derived from C under the gauge-fixing conditions 
N"" — 0. The answer, from a general perspective which 
is summarized in Appendix can be found in jl^ , and 
in general is in the negative. 

Bringing the results of to our case, we find (here 
we use [C] for the Euler-Lagrange variation) : 

[£] =0 , AT'' = 0^ [/:gf] =0 , Tia =0 , (6.1) 



7 



Pons, Shepley 



Bianchi-Type Cosmology 



To appear in Phys.Rev. D15-3 



where Tia are the Hamiltonian momentum constraints 
in the anholonomic basis and expressed in coordinate- 
velocity space. The N°- have been set to zero in Ha- If 
we also implement the Killing conditions in Tia, we get 



N 



T^a = -H^{gbc\a - 9ac\b)9"' 



(6.2) 



where A:^ = g'"''gdc- 

Therefore, if we use the completely reduced La- 
grangian, with gauge fixing plus Killing conditions, with 
configuration variables N, gab, we must be aware that the 
correct equations of motion will require that initial con- 
ditions be taken such that 



^bc'^a 







(6.3) 



We show in Appendix |^ that if the initial conditions sat- 
isfy equation ( |6.3| ), then equation (6.3) will hold at any 
time, provided the equations of motion arc satisfied. 

In Class A models, where C^^ — 0, equation ( |6.3| ) has 
only the term C^i^k'^ = 0. It is worth noting that, except 
for Type I models, which have vanishing structure con- 
stants, in all other Class A models, equation (^^) must 
be enforced on the initial conditions. A diagonal form 
for the metri c, fo r instance, will guarantee the fulfillment 
of equation ( |6.3| ) for most Class A models as long as the 
structure constants are taken in the form displayed in p^ ] 
(the exception is the group of Type VI_i). However, this 
diagonal form may not exhaust the possible physics avail- 
able in these models, and we do not assume diagonality 
here. 

Notice that the HPD cann ot b e used to make the initial 
conditions satisfy equation ( |6.3| ): Under a change 



gab 



g'ab 



the momentum constraints change as 

Ha(g') = (detM)M,^Hfc(g) . 

The result of our analysis is the following: The com- 
pletely reduced Lagrangian gives the correct equations 
of motion for Class A models if and only if we tune the 
initial conditions (at t = 0) gab{0), 5afc(0) in such a way 
that 



C^.g'-'m.M = 



(6.4) 



To this result we must add the fact that two three- 
metrics, gab and g'ab, ^'i'^ physically equivalent, that is, 
gauge related, if there is an HPD such that in the no- 
tation of equation (5.1) g'ab = MagaiM^. We must im- 



plement this fact in counting the number of independent 
initial conditions. In order to have a correct dynam- 
ics and a correct counting of the degrees of freedom for 
Class A Bianchi models we must require that (1) t he i ni- 
tial conditions must be chosen such that equation ( |6.4[ ) is 
satisfied, and (2) initial conditions related by HPD must 
be considered physically equivalent. 



VII. HAMILTONIAN APPROACH 



The standard Hamiltonian ADM approach [|T5[ is built 
upon the Lagrangian equation (^^) . Since C does not de- 
pend on N and N°-, the conjugate variables Vo and Va 
to N and iV" are the primary constraints in phase space. 
The canonical Hamiltonian has the form 



/ d^x{NHo + N^Ha) 



(7.1) 



The Dirac Hamiltonian iJn is formed by adding to H 
a linear combination (with arbitrary functions Aq, Aa) of 
the primary constraints: 



Hr 



H + / d3x(A°7'o + X^Va) 



(7.2) 



Stability of the primary constraints under the evolution 
generated by Hd leads to constraints which are now sec- 
ondary: 



Ho^Q : Ha^O 



(7.3) 



and no more constraints arise. 

All eight constraints are first class, and the four in- 
dependent gauge transformations that account for the 
diffeomorphism invariance are made out of them. For 
the sake of completeness let us write the gauge genera- 
tors which act through the Poisson bracket in the whole 
phase space (three-space integrations are understood 
for all repeated indices): 



(7.4) 



where the functions D'^,„ describe the first class structure 



of the secondary constraints, 

{Tip, Tip} = D'l^pHi, 



(7.5) 



and ^'^ are arbitrary functions. The relationship of the 
transformations generated by G{t) to the standard dif- 
feomorphisms generated by a vector field e^c?^ is given 
by (see [H for instance): 



Hi j\Ta 



(7.6) 



On a given metric, every infinitesimal spacetime diffeo- 
morphism is matched by G{t) with the specific func- 
tions ^'^ dictated by equation ( |7.6| ). This means that 
we have exactly the same gauge group either in the La- 
grangian or in the Hamiltonian formalism, here extended 
to the phase space, and that in this last formalism the 
same steps must be taken to reduce the gauge group to 
the time-reparameterization invariance. 
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The implementation of the Kilhng conditions and the 
gauge fixing TV" = into the Hamiltonian equation ( |7.2| ) 
leads to the same conclusions we have arrived at in the 
Lagrangian formulation. In particular, only for Class A 
models is the implementation of the Killing conditions 
into the Hamiltonian correct, in the sense that it gives 
the correct equations of motion. The reason is that the 
functional derivatives in the Hamiltonian formalism de- 



velop terms like the ones in equation (3.(:). In fact, any 



whereas the reduction of the right side of equation ( |7.8| ) 
is, after factoring out the three-space volume, 

i^c njrcd 

It is only for C lass A models that we get this result in 
equation (7.10): Only for Class A models do the imple- 
mentation of the Killing conditions and the computation 
of the Poisson brackets commute. 



of the Ti,^ is of the form 



Tip = |a;|Hp , 

in such a way that the three-space functional deriva- 
tives, corresponding to the Hamiltonian equations of mo- 
tion, become (here X stands for gaf,,7r''^, where is 
the canonical conjugate of gab in the reduced formalism: 



5X 



^^--(Y. + C^J(^) 



dX 



'dYaX' 



a>b 



dYaYbX 



(7.7) 



(As a matter of fact, only Hq depends on second order 
space derivatives.) 

To illustrate this last result, let us write the Poisson 
brackets for the momenta Ha- They satisfy 



(fxHaQ.jd^xHbi, 



(7.8) 



When we implement the Killing conditions on the quan- 
tity / d'^xHaS."', in order to get rid of the three-space 
dependence, it seems that we are bound to take the func- 
tions as constants (but this is not exactly true, as we 
will see shortly). This would generate diffeomorphisms 
of the type equation (^^) , namely inner automorphisms 
for the left-invariant Lie algebra. Then Tia reduces to 



n 



red 



^ab^ 9dc 



Cb ;;^cd 
be 



9da 



(7.9) 



(this is the phase space version of equation ( ^.2| )). 

To continue, we must first introduce the Poisson brack- 
ets for the reduced theory. The reduced Poisson brack- 
ets {— are defined through a renormalization of 
the old Poisson brackets {— ,— = F{— ,— }, where V 
is the (perhaps infinite) right invariant volume element 
V — J u}^ A Lj'^ A iv^ = J d^x\uj\. Notice that as long 
as the Killing conditions hold we can invert the relation 
tt"'' = \uj\n'''' to n"'' = j^J d^xir"^. The commutation 
relations for TVi?'^ are 



i'^a > /^b \r- 

'^^dc^^ab^'^^ 



red 



9fe 



'9eb 



r<f ~c 



(7.10) 



A. Homogeneity-preserving diffeomorphisms, 
revisited 



We have said, after equation ( [7.8[ ), that the imple- 
mentation of the Killing conditions on the generator 
/ d'^xHa^'^ does not require the functions to be con- 
stant. Th e re ason is as follows. In the computations in 
equation ( |7.8| ) we have been dropping boundary terms. 
This is correct if our functions 5" have a compact sup- 
port on Ej. But we know from the Lagrangian analy- 
sis that some diffeomorphisms that do not vanish at the 
boundaries, for instance the HPD, play an important role 
when we perform the dimensional reduction through the 
Killing conditions. 

Consider the generator of three-space diffeomorphisms 
for the reduced set of variables gab,'^"'^, 



I 



d^xHaC 



with 



In order to get the standard result 



Sgab 



gab, J d^xT-CaC^ ^^b\a+L\b 



(7.11) 



(7.12) 



(7.13) 



we a lways perform integration by parts in equation 
( 7.11 ). Since the generator equation ( 7.11 ) acts locally, 
what we actually do in practice is to compute Sgab{t,x) 
using functions that are identical to ^a in a neighborhood 
of a given point but that vanish at spatial infinity. If 
we pe rform directly the integration by parts in equation 
(7.11) we get 



(7.14) 



which differs from equation (7.11) at most by a boundary 
term. 

Boundary terms may or may not exist depending on 
the topology of E. As long as gauge tra nsfor mations on 
local functio ns ar c considered, equation (7.14) is as good 
as equation (7.11) to generate the transformations equa- 
tion (7.13), but Gs is better when we must work with 
arbitrary functions £,a that are different from zero at the 
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boundaries. Notice that whereas equation (7.11) van- 
ishes on a solution of the eq uation s of motion, because 
Ha is a constraint, equation (7.14) does not need to do 
so in the case when the functions do not vanish at the 
boundaries. Thus G3 is not bound to be zero on a solu- 
tion of the equations of motion. It is not a constraint, 
except when there are no boundaries in E or for some 
other some particular cases that will be examined below. 

The dimensional reductions of equation ( 7.1l| ) and 
equation ( [7.14 ) give different result, because integration 
by parts for the space coordinates does not exist in the 
reduced formalism. Now we apply the Killing conditions 
to equation ( [7.14 ): 



where the connection coefficients F^^ are 

(here gf'^ = e^"^ since N" ^0). Therefore, 

G3 - 1 d^x\u;\n'^gUYb(^ " ■ (7-15) 

Dimensional reduction requires us to factor out the vol- 
ume V = J (Px\uj\. This means that we need the remain- 
ing piece in the integrand, TT:'"^gcd0^bi'^ — C^gC^), to be 
independent of the space coordinates. Since tt'"^ and gcd 
are already space independent because they satisfy the 
Killing conditions, we end up with the requirement 



Bi 



(7.16) 



with a constant matrix. This condition is exactly 



equation (5.2), the condition for HPD, expressed in the 
dual basis. Thus we see that the dimensional reduction 
in phase space has room for and only for the HPD we 
already found in the Lagrangian formalism. Then G3 be- 
comes G3 = VB^TT^'^gcd, and since the reduced Poisson 
brackets are = ^{^i we get for the canoni- 

cal generator of HPD in the reduced formalism. 



G3 



9cd 



(7.17) 



We find here another special feature of the Bianchi 
models, for the G3 in equation ( [7.17 ) generates gauge 
transformations, HPD, and yet it is not a constraint of 
the reduced formalism by itself, except for the case in 
Class A models of inner automorphisms. In this last 



case, we have B 



d 



B'^C^, and G3 for Class A mod- 



els becomes G3 = B°'W^^, with W^'^ given in equation 
( [7.9| ) with G^g = 0. These results completely match the 
ones previously obtained in the Lagrangian analysis. 

The generators of HPD that define automorphisms of 
the Lie algebra that are not inner automorphisms (called 
"outer" HPD in |^) are always constants of motion j|] of 



the reduced formalism. They are enforced to be con- 
straints if S has no boundaries, although these con- 
straints cannot be retrieved as such constraints from the 
reduced formalism by itself. In any case, either being 
constants of motion or constraints, they always generate 
gauge transformations in the reduced formalism. These 
transformations, in either the Lagrangian or the Hamilto- 
nian picture, are available as gauge transformations from 
the outset, as a consequence of the rationale of the re- 
duction procedure, but they cannot be retrieved as gauge 
transformations from the reduced formalism by itself. 
A comment is in order for Class B models. The gen- 



erators equation (7.17) give the HPD for Class B as well 
as Class A models. In particular, C^f^n^" gcd generates 
the inner automorphisms associated, in the configuration 
space picture, with the right-invariant vector Ya . Never- 
theless, G^^7r'"^gcd is not necessarily a constraint for the 
Class B models because the true momentum constraints 



are those of equation ( 7.9), as one can directly verify from 
the equations of motion. Oddly enough, in Class B mod- 
els at least one of the generators of the residual gauge 
transformations associated with inner automorphisms is 
not a constraint. Also, according to equation (7.10), 
these constraints, two of them, are no longer first class. 
These results hint that the Lagrangian or Hamiltonian 
formulation for Class B models is not possible in general. 

Notice that in the case of three-manifolds S with no 
boundaries, the possible new constraints G3 in equation 
(7.17) appear also in the Lagrangian formalism by simply 
writi ng do wn the process from equation (7.11) to equa- 
tion (7.14) in configuration- velocity variables. 

One could think that the mechanis m of integration by 
parts , as used to transform equation ( 7.11 ) into equation 
(7.14), and that helped us find the whole set of HPD in 
phase space, could perhaps be used in Class B models to 
transform the reduced Hq into the correct generator of 
time translations. In fact, some unwanted pieces can be 
eliminated this way, but this does not solve the problem. 



B. Degrees of freedom 

Notice also that the tangency of the Hamiltonian evo- 
lution operator equation (7.1) to the gauge fixing sur- 
face defined by N"' — implies A° = 0. On the other 
hand, once the gauge fixing N'^ = has been intro- 
duced into i?D, the stability of the primary constraints 
Va = will no longer produce the secondary constr aint s 
Ha — 0. This result parallels that of equation ([6j|). 
As a matter of fact, and since they have become second 
class constraints, these three couples of canonical vari- 
ables iV°,'Pb, are readily eliminated from the formalism 
by taking the Dirac bracket, which is nothing but the or- 
dinary Poisson bracket for the rest of the variables (the 
same number of variables, N'^, are eliminated at the 
same stage in the Lagrangian formulation). The impor- 
tant point is that the momentum constraints are gone 
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from the reduced formalism, but yet they must be imple- 
mented (in case they do not vanish identically) from the 
outset as restrictions imposed on the initial conditions, if 
we want to have the right equations of motion. 

It is remarkable also that the generators of HPD in 
the dimensionally reduced theory are not necessarily con- 
straints, except for the case of inner automorphisms and 
Class A models. We still need to fix all the residual gauge 
freedom corresponding to the three-diffeomorphisms, be- 
cause these three-diffeomorphisms have been with us 
since the beginning. We are in the same situation as we 
were in the Lagrangian formulation, and therefore the 
same HPD, exactly the same, appear here, as we have 
just shown. That is why the degrees of freedom in both 
the Lagrangian and the Hamiltonian formalisms coincide. 

In 1^ a discrepancy between Lagrangian and Hamil- 
tonian degrees of freedom is argued on the basis that 
the only HPD available in Hamiltonian formalism are 
those yielding inner automorphisms. Our analysis dif- 
fers from the one in in that we show that the resid- 
ual three-diffeomorphism invariance that still needs to be 
fixed corresponds to the HPD that we found earlier in the 
Lagrangian formalism. And this is something that one 
knows in advance through the process of the reducing the 
original gauge group. 

The claims in |^ are opposed in |^ , where it is argued 
that the degrees of freedom for both the Lagrangian and 
Hamiltonian formalisms are the same. However, they 
base their claim on rejecting the HPD that are not asso- 
ciated with inner automorphisms as residual gauge trans- 
formations that need to be fixed. In the next section we 
show that all the HPD are indeed residual gauge trans- 
formations that describe redundancy in the initial condi- 
tions of the system and therefore do need to be fixed. 



VIII. CONCLUSIONS 



tion of the gauge fixing that sets the shift variables to 
zero. As we have just said, the first problem prevents 
the Class B Bianchi models from having a reduced La- 
grangian or Hamiltonian formulation in general. The 
second tells us that the initial conditions must be cho- 
sen to satisfy some former constraints (the momentum 
constraints) of the original theory, even though they are 
no longer constraints for the reduced one. 

We have shown that the Hamiltonian formalism has 
the same gauge freedom that is available in the La- 
grangian formulation. The reduction of the gauge group 
in phase space follows exactly the same steps as in 
configuration-velocity space, and the same features, the 
same problems, the same considerations, and the same 
results apply as well. The fact that the reduced momenta 
Ha in Class A models can only generate HPD that are 
inner automorphisms of the Lie algebra has its counter- 
part in configuration-velocity space in that these HPD 
are the only ones that can be derived as gauge transfor- 
mations from the reduced Lagrangian on its own. But 
that does not mean that the HPD associated with the 
"outer" automorphisms, if there are any, do not need to 
be quotiented out. Rather, all HPD must be quotiented 
out, either in the Hamiltonian or in the Lagrangian for- 
malism. 

Also, what is and what is not a gauge symmetry and 
the role of the residual gauge symmetries has been ana- 
lyzed in detail. We think that our considerations throw 
a definitive light on some issues that have not yet been 
settled in the literature p|j20[|. 

Through our analysis, the gauge group is reduced in 
several stages, and residual gauge symmetries appear. 
According to the fact that they all proceed from the 
original gauge group, we claim that these residual gauge 
symmetries must be fixed thoroughly. We think that this 
is an important point that deserves further elaboration 
below. 



We have proved that Class A Bianchi models allow for 
a Lagrangian formulation as a mechanical particle-like 
system (finite number of degrees of freedom, in contrast 
to field theory) as long as the setting of the initial con- 
ditions is taken in a certain way. 

In our transit toward this result we have proceeded as 
follows: We have carefully produced a partial gauge fixing 
in four steps, to reduce the initial gauge group generated 
by the four-diffeomorphism invariance to time reparam- 
eterization invariance. As a consequence, the troubles 
that beset Class B Bianchi models are identified as ob- 
structions to the commutativity between two processes: 
implementing the Killing conditions either into the equa- 
tions of motion or directly into the Lagrangian, followed 
by deriving the equations of motion. 

We distinguish, as regards the Lagrangian or Hamil- 
tonian formulations, two types of problems, which corre- 
spond, respectively, to the implementation of the Killing 
conditions into the Lagrangian and to the implementa- 



A. Gauge fixing for tiie residual gauge group 

Let us show the need to fix the gauge for the resid- 
ual gauge group of three-space diffeomorphisms that ap- 
pears at the end of the third step of the gauge fixing 
procedure (Section IV). Consider that we have a solu- 
tion gij{x, t), N{t) of the Einstein equations with A^' = 
and with initial conditions at t = 0. Now consider the 
action of the infinitesimal diffeomorphism transformation 
defined by the vector field e'^di, 

e\x,t)^ fit)i\x), 

where / has been chosen to satisfy 

/(O) = , /(O) = ; fit) = 1 , fit) = for < > 1 . 

Since this vector field generates a foliation-preserving 
diffeomorphism, it will define a gauge transformation 
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which is still allowed in the formalism after the first- 
step gauge-fixing has been done. The transformed met- 
ric g'ijix), N'{x), N'\x) is also a solution of the Einstein 
equations. Notice that both metrics and 5'^^ share 
the same initial conditions at i = 0. It is then obvious 
that the correct interpretation of a gauge transformation 
dictates that g^^ and g' are physically equivalent and 
must be so for any other time we take for the setting of 
the initial conditions. 

Now consider the relation between these two metrics 
but with initial conditions taken at i = 1. They are re- 
lated by an infinitesimal three-space diffeomorphism gen- 
erated by the time-independent vector field e^{x)di. The 
two metrics satisfy the gauge N''' = 0. The gauge trans- 
formation that connects them is not localizable in time 
(being different from the identity only within a finite time 
interval), which is the point made in pO[ | to deem this 
transformation as non-gauge, and yet both metrics must 
be physically identified. This result proves that we must 
fix the gauge even for the three-space diffeomorphisms, 
namely, the residual gauge group. It is in this place that 
the HPD play the significant role we have seen in Section 

and Section [vi^ . 

Our analysis agrees with that of |^ with regard to the 
Lagrangian formulation, and we agree with their space- 
time counting of the degrees of freedom (see the table in 
1^ ; also see p^-[l8|| for a counting of degrees of freedom in 
models with compactificd spatial sections) . But we differ 
in other respects. Let us make the differences clear; they 
concern the status of the HPD associated with "outer" 
automorphisms of the right invariant Lie algebra. 

In g, "outer" HPD, that is, HPD yielding automor- 
phisms of the Lie algebra that are not inner automor- 
phisms, count as gauge degrees of freedom in the La- 
grangian formulation, but are considered non-gauge sym- 
metries in phase space. In these "outer" HPD are 
always taken in both formulations as non-gauge symme- 
tries. Instead, from our systematic procedure of reducing 
the gauge group, we deduce that the "outer" HPD are 
gauge transformations in both the Lagrangian and the 
Hamiltonian formalisms, and they always count as gauge 
degrees of freedom. 

Also see our comments in Appendix For exam- 
ple, we explicitly consider the case of the Bianchi Type 

1 model, where the number of degrees of freedom is 1 (if 
the surfaces of homogeneity have the topology of R'^). It 
may appear to some people that an odd number of de- 
grees of freedom in a Hamiltonian formalism is somehow 
not correct, but that is the case here, and we discuss this 
matter in a bit more detail in Appendix ^ 

Summing up, we have exhibited in a simple way the 
two different problems that appear when we implement 
the Killing conditions and the gauge fixing into the La- 
grangian or into the Hamiltonian. The first problem 
prevents the Class B Bianchi models from having La- 
grangian or Hamiltonian formulations (except possibly 
in special cases), whereas the second is solved by intro- 



ducing some requirements on the initial conditions. We 
have also shown that there is no ambiguity or matter of 
interpretation in what must be understood as a gauge 
transformation, either in the Lagrangian or the Hamil- 
tonian formalism. Finally our analysis proves that the 
number of degrees of freedom in both formalisms is al- 
ways the same. 
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APPENDIX A: GENERAL SPATIALLY 
HOMOGENEOUS METRIC 

In order to demonstrate some of our ideas explicitly, 
we here give formulas for the connection, curvature, and 
field equations for a general spatially homogeneous met- 
ric, including a general lapse function and shift vector. 
A word on calculational procedure: It is easiest to use an 
orthonormal basis adapted to the Killing structure: 



where the basis is defined by 

(T° = Ndt , cr* ^ 6j,(a;° + N^dt) 



(Al) 



(A2) 



where iV, A^", gab are functions only of t, and we use the 
convention that bold face denotes a tensor or a form; the 
one-forms a;" obey the following relation, indicating that 
they are invariant under the isometry group whose lie 
algebra is defined by the structure constants C^^: 



Note that the three-metric is defined by 



gab 



(A3) 



(A4) 



and we will consistently use for orthonormal in- 

dices, and a,b,c, . . . for ordinary three-indices. The ma- 
trix B — (bait)) represents an arbitrary square root of 
the metric; precisely what form it takes will be irrele- 
vant. The inverse of the three- metric is defined by 



ac ca 

e gcb = Ofc 



(A5) 



in order not to confuse it with the ah components of the 
contravariant four-metric. The inverse of B is the matrix 
A = (af), so that 



where a%l, = 



atbl = 5t . (A6) 
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The first Cartan equations witli tlie torsion set equal 
to zero are 



The Ricci tensor components, the scalar curvature, and 
the Einstein tensor components are defined by 



(A7) 



where cr(^ are the connection one-forms, used in forming 
the covariant derivative of a tensor. The metric compati- 
bility equations (covariant derivative of the metric equals 
zero), when the metric components are constants, as in 
an orthonormal basis, here read 



(A8) 



Equations (A7) and (A8) uniquely determine the connec- 
tion one- forms. The second Cartan equations are 



(A9) 



they determine the Riemann tensor components (see p5| ] 
for more details). 

It is convenient to define a matrix related to the loga- 
rithmic derivative of B and its symmetric and antisym- 
metric parts as the matrices K, L, M. It is also conve- 
nient to define the orthonormal projection of the struc- 
ture constants (which then become a time-dependent ar- 
ray D plus another array E by (here ' denotes d/dt): 



The evolution equations (in vacuum) are setting the 
space-space components of the Ricci tensor to zero, where 



Rj — J^^^^ ~ LikMjk — LjkMik + LkkLij 



-{D],+Dl,)D 



(A15) 



Next we display the time-time and time-space compo- 
nents of the Einstein tensor; these set to zero are the 
constraint equations on initial value data: 

Gq — -LstLst — -^{Lss)'^ 

+ IdI,D^^ + il?L^st + ^^L^L , (Al6a) 



(A16b) 



(AlOa) 



^ij — 2^^^i ~^ ^1^^ ' ^^^^ 



-(K,j-Kj,), (AlOb) 



Finally, we will need the scalar curvature in order to 
display the reduced Lagrangian. In the action integral 
the volume element is really 

(T° A cr^ A A (T^ = iV^di Ao;^ Aw^ Aw^ , (A17) 



D)k = K^H^Sc : E^jk - D), -Dl^-D%. (All) 

With these definitions, we have for the curls of the or- 
thonormal basis forms 

dcrO = , dcr' = K.jCT^ A cr^ + ^OLa^ A . (A12) 

2 

The connections forms are: 

<tO ^al^ L,^a^ , a] = -M^.tr" + ^E.^k^'^ . (A13) 

The results for the independent components of the Rie- 
mann tensor are most conveniently displayed after raising 
an index (using S^^): 

Rq] — J^^-^ij + LikLjk ~ LikMjk ~ LjkMik , (Al4a) 



Rm — EijD^^ + -{LjkEj.i - LjiEjik) , 



2 

iiLjk 

1 „ „ 1 „ „ 1 



(A14b) 



R^i — LikLji — LifL 



^EimkEjml + -^EimlEjmk + T^EijmD^^ . (A14c) 



where g = Aei{gab), and the spatial integral can be set to 
the constant V: 



The action integral is therefore 



X= / LVdt 



where 



= {2^Lssr-2^N^hlDl,L,s 
+N^(^LstLst-{Lssf 



(A18) 



(A19) 



and note that we have separated out a total time deriva- 
tive. Here it is best to be explicit, and we display the 
Lagrangian in terms of the configuration space variables 
N,N-,gab: 
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C 



N 



+ lg,,e^''N'^Nfc:,C% + lN'^Nfc:,Cl 



-N^g 



gadS e Of,j,Og^ 



(A20) 



The vacuum field equations are most conveniently 
written as four constraint equations and six evolution 
equations: 



G[] = , G° - , i?; = 



(A21) 



These equations can be derived from the Einstein-Hilbert 
Lagrangian density before imposing any symmetry re- 
quirements and then imposing the symmetry require- 
ments after the equations have been derived. We will 
compare these equations with those derived from the re- 
duced Lagrangian £; in other words, by imposing the 
symmetry requirements first and then deriving equations. 

The equation 6C/SN = is readily seen to be exactly 
Gq = (up to a factor of 2y/g). 

The equations S£/6N°- = (after canceling N and y^) 



are 



= +gdbg'''Ct + 2Cl,N'Ct 



bd 



+ geig-'N'cifil, -f N'^CtCl, 



(A22) 



This is not proportional to the G\ equation in a Class B 
model; the difference (up to nonzero factors) is 



{gabg"^ + 2G,^,iV/ - ga.g^^N'^CU C^^ 



(A23) 



In other words, if the lapse is kept in, then all of the 
constraint equations can be derived for the Class A mod- 
els. In Class B, if these terms happen to be zero any- 
way, well and good, but they won't automatically van- 
ish. Note several things about the terms in equation 
( [A2^ ): The first one is independent of 7V°; therefore, one 



might hope that adding an appropriate term to the La- 
grangian might remove it (a term, that is, linear in iV"). 
The second term has N"" multiplied by a matrix which 
is symmetric; supposedly it too could be removed with 
an appropriate term in the Lagrangian. The last term, 
however, has -/V" times a matrix which is not symmetric; 
it would seem that there is no hope of generating this 
term by adding something to the Lagrangian. 



APPENDIX B: SPATIALLY HOMOGENEOUS 
MAXWELL THEORY 

A relatively simple system which illustrates many of 
our ideas is that of a spatially homogeneous electromag- 
netic potential which obeys the Maxwell equations. (Here 



we start with the four- vector potential rather than with a 
spatially homogeneous field tensor, which of course could 
be produced using a non- homogeneous potential vector.) 
The background metric is taken to be a simple spatially 
homogeneous one, so we do not require that it obey any 
particular field equations. We take the metric to have 
components 77^,^ in an invariant basis: 



where 



2 ^ 



(Bl) 



(B2) 



Note that there is a freedom to transform the spatial 
basis {w'} by an orthogonal transformation but not nec- 
essarily by a general linear transformation. 

The vector potential actually is an equivalence class of 
one-forms related by gauge transformations. Here we de- 
mand that the class contain at least one member which 
is spatially homogeneous. If that one is expressed in the 
invariant basis, its components are functions only of t: 



A = An{t)dt + Ai{t)ijj^ , 
and it produces the following field 2-form: 



F = dA = A,dt A 



iA,Gj,a;^- Ao;^ 



(B3) 



(B4) 



Note that Aq has disappeared; it plays no part either in 
the Maxwell equations or in the Lagrangian. It can be 
made to vanish by a gauge transformation of the kind 
A — * A' = A + dA(i). 

We choose to concentrate on the metric as expressed 
in an invariant basis; this is a choice of metric gauge. 
Maxwell theory cannot be expressed without at least 
some reference to a background spacetime metric, and 
so this metric gauge affects the electromagnetic poten- 
tial. We also choose to concentrate on the vector po- 
tential within its equivalence class which explicitly has 
vanishing Lie derivative with respect to the generators 
of the invariance group. This, too, is a choice of gauge, 
though as we discuss below, in some cases there is still 
some residual gauge freedom. 

The easiest way to calculate the Maxwell equations (in 
vacuum) is first to write down the dual field 2-form: 



"F = ie.jfe^G^fcdt Aw* - ieysi^a;* A uj^ , (B5) 



where e^fe is the Levi-Civita symbol, equal to ±1 if {ijk) 
is an even/odd permutation of (123) and to zero other- 
wise. The Maxwell equations are: 



= d*F 
1 



tstuAs + i^ezjkAsC]kCl^ dt A A a; 



-eiskAkCl^u}' A u;' A u}" 



(B6) 



14 



Pons, Shepley 



Bianchi-Type Cosmology 



To appear in Phys.Rev. D15-3 



and they are equivalently written as: 

, (B7a) 

= + lAsq.q, + A.CttCr^ . (B7b) 

One can double-check these equations by forming the 
connection one-forms and computing F^";„ — 0. 

The action integral for a general Maxwell vector po- 
tential is is the determinant of the spacetime metric 
in a coordinate system): 



(B8) 



In this case, VtsI independent of t, and so the reduced 
Lagrangian Cr is: 



1 



+ i^^q^C], . (B9) 



The Euler-Lagrange equations = 6Cii/5Ai are: 
= + -AsCjf.Cji. . 



It i s clear that in Class B models equations (B7) and 
( ^10| ) do not agree. There are two differences: Equa- 
tion (B7a) is a constraint equation, and it simply cannot 
come from a variational principle which is homogeneous 
quadratic in the velocities. Equation ( B7b| ) has the term 



AsC-^Cf^ which do es no t appear in equation ( |B1C| ). No- 



tice that equation (B7a) is automatically satis fied i n the 
Class A case (CJ^ = 0), and also equations (B7b) and 
( plO| ) then do agree. 

One thing wh ich is true, however, is that the con- 
straint equation ( B7a) is compatible both with equation 
(B7^) and equation ( B10|). To see this fact, take the time 
derivative of equation 07a ): 



— AjCjs 



a homogeneous one-form whose curl vanishes (and there- 
fore which can at least locally be expressed as the curl of 
a function). Let this one-form be 



(Bll) 



where the components k^j are functions only of t. We 
require that the curl of k be zero: 

= d/« = k^dt A -f ]-K,CL^^ A US'" . (B12) 
2 

Therefore kq is arbitrary, and ki must be a set of con- 
stants subject to the condition 



(B13) 



Only Bianchi Types VIII and IX require that be zero. 
In fact, in a Class B model, Ki may be taken to be pro- 
portional to Cfj . 

We illustrate with the example of the general Class 
B model, in which only Cf^ 7^ (this prescription can 
always be satisfied in any Class B model by using an 
orthogonal transformation of the invariant basis). The 
Jacobi identity then requires that Clj — 0. The reduced 
(BIO) Lagrangian in this case is 

Cr = -i(i? + AI + Aj) + \aaAbC%C% , (B14) 

where the indices A, B range only over (2,3). The Euler- 
Lagrange equations are 

1 



where 



^ , = + -^AbMab 



Mab = Cfi^Cf,, 



When equation (B7b) is multiplied by C*,, the result is 



The Maxwell Equations ( |B7| ) read 

= ii , = + \AbMab + AbNab , 



— C^isA-i + -^AtC^^^Cis + AjClfi^^Clg . 

The middle term vanishes because of the Jacobi identity, 
which in a three-dimensional Lie algebra is equivalent 
to = Cf C*(. The last term, which would be absent 



where 



anyway if equation (BIO) had been used, vanishes as a 
consequence of the antisymmetry of Cl^ in its lower in- 
dices. 



Thus, even if the constraint equation ( B7a ) were put in 
by hand, the problem would remain whether any varia- 
tional principle could reproduce the last term in equation 



(ill)- 



We now turn to gauge transformations. A gauge trans- 
formation here is the addition to the vector potential of 



Nab = CaiCIs 



The reduced Maxwell equations thus differ in two ways 
from the Euler-Lagrange equations: First, the Maxwell 
equations have the constraint equation Ai ==0. The 
residual gauge freedom allowed by equation (B13) says 
that Ai can be made zero by a choice of gauge. In con- 
trast, the Euler-Lagrange equations only require that Ai 
be at most a linear function of and the gauge freedom 
only will allow the value of Ai to be set to zero at a 
particular time. It is not, however, difficult to add the 
constraint = to the Euler-Lagrange equations in an 
ad hoc manner. 



15 



Pons, Shepley 



Bianchi-Type Cosmology 



To appear in Phys.Rev. D15-3 



Second, the Maxwell equations for Aa have the addi- 
tional term involving Nab- If it happens that Nab is 
symmetric, then a Lagrangian can be found to repro- 
duce this term: The Lagrangian would add the term 
^AaAbNab to the reduced Lagrangian equation (B14). 
However, if Nab is not symmetric, then such a La- 
grangian would not in general be possible — certainly 
no Lagrangian could be found to reproduce the Maxwell 
equations exactly, though under some circumstances it 
may be possible to find a Lagrangian which would pro- 
duce equations equivalent to the Maxwell equations. 

For example, the standard structure constants of a 
Bianchi Type III model |l| are = -C21 = 1^ the 
other structure constants vanishing. In this case Nab = 
diag(l,0), and the Maxwell equations can indeed be de- 
rived from a variational principle (with the constraint 
= Ai being put in by hand): The Maxwell equations 
are: 

= Ai = Ai = A2 ^ A3 , 

which clearly can come from a constrained variational 
principle. Note that in this Type III case the reduced 
Lagrangian yields as Euler-Lagrange equations: 

= ii = + ^2 = ^3 , 

which are not at all the same as the Maxwell equations. 

A second example is that of a group which has struc- 
ture constants = —C21 = Cf^ = — Cf^ = 1, the 
rest being zero. In fact, this group is also Type III, in 
a basis which is a linear transformation of the basis in 
the preceding example. Since this transformation is not 
an orthogonal one, the Maxwell equations in this case 
differ significantly from the preceding. The reduced La- 
grangian is 



Al 



The Euler-Lagrange equations are: 

= ii = + 2^2 = ia ■ 
The Maxwell equations are: 

Q = Ai= Ai= A2= A3- A2 . 

The difference between the Maxwell equations and the 
Euler-Lagrange equations in this example are profound: 
First, the Maxwell equations include the constraint Ai — 
0, as in the previous example. Second, the Maxwell equa- 
tions cannot be derived from a variational principle, un- 
like the previous example. 



APPENDIX C: RELATIONSHIP BETWEEN 
SETTING THE GAUGE AND VARIATIONAL 
PRINCIPLES 

Here we summarize and expand the main result of p4[ . 
For simplicity we will use the language of mechanics (fi- 
nite number of degrees of freedom), although everything 



can be translated to field theory. Consider a singular 
Lagrangian L that leads to some primary constraints, 
<j)p — Q in phase space, presumed to be effective (each 
has non- vanishing gradient on the constraint surface). 
For the sake of simplicity we will consider the case when 
all the constraints, primary, secondary, etc., are first class 
(this is what happens in our generally covariant theories). 
The equations of motion obtained from L are: 



5L_ 

5q^ 



where 



and ai 



tjq 



dL 



■J — 



dq'-dq^ 



Wdt 



The Lagrangian dynamics can equivalently be de- 
scribed by a vector field that exists on, and is tangent 
to, the constraint surface in configuration- velocity space: 



X 



d_ 

di 



dq'- 



a\q,q) 



d_ 



(CI) 



The a* are determined from the equations of motion and 
the stabilization algorithm; rjP are arbitrary functions of 
time (or spacetime in field theory) and any other variable; 



and Fp is 



:_d_ 

'dq' 



where 7^ are the null vectors of Wij . These null vectors 
can be given as 



dp., 



(C2) 



where Pi( q, q) = dL/dcf'. Let us point out that ?yTp in 



equation (|C1) is the piece that corresponds in our case 
to equation (lT|). 

Notice that 0^7^ = is a consequence of the equa- 
tions of motion SL/Sq' — 0. They are called the primary 
Lagrangian constraints. 

Consider now a partial gauge-fixing of the dynamics 
given by a set of new constraints Xp' ~ 0, with \p'\ < 
defined in configuration space (holonomic constraints). 
Let us split the set of indices p into two sets, p' and p", 
in such a way that Tpux^i = and iFp/Xo-'l 7^ 0. Then, 
the requirement X(xp') = determines the functions 77'' 
and leaves completely undetermined rjP . 

In such a situation, the following result is proved in 
[Q: If we plug the holonomic gauge-fixing Xp' = into 
the original Lagrangian L to get the reduced Lagrangian 
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Lqp, then the following equivalence holds (where [L] is 
the Euler-Lagrange variation): 

[L]=0, xp'=0^ [Lgf] = , a,7;, = , (C3) 

where ctiYp, are a subset of the primary Lagrangian con- 
straints for L with the understanding that the gauge fix- 
ing has been plugged into them. 

So we see that to describe the same motions in the re- 
duced space it is not enough in general to impose the new 
Lagrangian equations of motion, [Lgf] = 0, but some ad- 
ditional constraints must be required, too. 

Now, we will prove, using the machinery of |l4[ , that 
in order for the solutions of [Lgf] = to satisfy the con- 
straints ai"fp, — 0, we only need to impose them in the 
initial conditions. 

Since the gauge-fixing constraints Xp' — are holo- 
nomic, they just reduce the configuration space. We have 
adopted the notation for the local coordinates in the 
original configuration space. We will adopt the notation 
for the local coordinates in the reduced configuration 
space. 

It is proved in [|lj] that the Lagrangian evolution op- 
erator in the reduced velocity space takes the form 



= XRo + ^''"fp" ; 



d 



(C4) 



where fjP are arbitrary functions; Tpn is 

d 



p" 



and similarly to equation ( |C2| ), 7°,, is defined by 

-(Q,P) , (C5) 



1> 



d(j)pn 



dPa 



where (fipn are the primary Hamiltonian constraints cor- 
responding to the reduced theory. It turns out [^ that 
these constraints are related to the original constraints 
(j)pn by 



(kp-'iliQ),?) = (f>p"{Q,P-g^) 



(C6) 



From all these results, the following equalities hold: 
fp"((a.7p')lq(Q)) 



: {Tp„ (a,7;,))|g(Q) 

: TL*{(l)p,,(l)p,>}\g(_Q) = , 



(C7) 



Here J-L* is the puUback of the Legendre map TL from 
velocity space to phase space. We have used in the last 



equality the fact that the puUback of a primary Hamil- 
tonian constraint is identically zero. 

With this new result, Tp" ((ai7p')|q(Q)) = 0, we can see 
that the arbitrary part in t he reduced Lagrangian evolu- 



tion operator equation (C4) has no effect on {ai"fpi)\q{Q)- 
This means that in the constraint surface for the reduced 
theory (this surface is determined by the equations of mo- 
tion for Lqp), the operator Xrq must also be tangent to 
the surface ((ai7p')lij(Q)) = 0- (This is the only way to 



ensure the equivalence in equation (C3), for if the opera- 
tor Xro is not tangent to the surface {{<^nl')\q{Q)) = 0, 
then there will be no solutions of the equations of mo- 
tion for the original Lagrangian in the gauge xp' — 0; 
and we know that these solutions exist.) Therefore, if 
we consider a solution of [Lgf] = with initial condi- 
tions satisfying (Q;i7*/)|q(Q) = 0, then the whole solution 
satisfies (Q;i7^,)|,(Q) = 0. 



APPENDIX D: DEGREES OF FREEDOM 

In this Appendix we illustrate with a very simple case 
the counting of degrees of freedom for our models. We 
also discuss the fact that spacetime rigid symmetries may 
be understood as residual gauge transformations. 

We start by counting the true degrees of freedom for 
the Bianchi Type I case. The counting can be done 
in configuration- velocity space or in phase space — the 
counting is the same [|l^. We start with the variables N 
and TV, N"- and 7V°, gab and gat- 20 apparent degrees of 
freedom. The shift vector variables may be eliminated by 
our general considerations, so that 14 apparent degrees 
of freedom are left. 

There is one constraint in the Lagrangian formalism. 
To fix the time reparameterization invariance we define 
the time parameter as a function of our dynamical vari- 
ables, excluding N. Its stability, that is, the fact that the 
time derivative of this definition must vanish, will give a 
new constraint, containing the variable N, with no time 
dependence. The stability of this last constraint gives a 
new constraint that can be used to isolate TV. The re- 
quirement of stability now determines the arbitrary func- 
tion in the Lagrangian evolution operator, and no more 
constraints appear. We are left with 1-f 3 constraints (one 
true constraint plus three gauge fixing constraints), low- 
ering the number of degrees of freedom to 10. The same 
counting can be done in phase space: In this case, there 
are 2 true constraints, and we must introduce 2 gauge 
fixing constraints. 

In Type I the momentum constraints are identically 
zero, so they do not reduce the degrees of freedom. We 
are only left with HP D sym metries. Every constant ma- 
trix of equation (7.16) defines an automorphism of 
the Lie algebra. 

Now topology enters the picture. If the topology of 
E is "open" (for example, M'^), simply connected, and 
with global Killing vectors, then to each matrix there 
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corresponds an HPD through equation ( |5.2| ). Therefore, 
there are nine HPD gauge degr ees of freedom. The gen- 
erators of HPD, equation ( [7.1'^), a re not const raints, be- 
cause in this case equations ( 7.11 ) and ( 7.14 ) differ by 
boundary terms. So the final number of degrees of free- 
dom for a Bianchi Type I model with surfaces of homo- 
geneity with the topology of M'^ is 10 — 9 = 1. 

It may appear strange to some people that a Hamilto- 
nian formalism can turn out to have an odd number of 
degrees of freedom, as in the above example. In fact, we 
have already said that HPD need only be implemented in 
the initial conditions, so in this sense we have a Hamil- 
tonian formulation plus an equivalence relation coming 
from the outset: It is not generated by the reduced for- 
malism but is a remnant of the generally covariant theory 
we started with. To find the number of degrees of free- 
dom number is a matter of counting constraints, which 
trajectories count as the "same" as others, and whether a 
parameterized trajectory or its orbit (its one-dimensional 
set of points in phase space) should be considered as 
physically significant. We have been guided by the prin- 
ciple that the Lagrangian and Hamiltonian formalisms 
should be equivalent (see ||l^) in coming to the conclu- 
sion that in fact they are. 

If the topology of S is that of a three-torus, then there 
are no infinitesimal HPD available (there are finite HPD 
but they do not reduce local degrees of freedom). We 
end up with a final number of 10 degrees of freedom. We 
will not consider other topologies. 

Let us turn to the second consideration in this Ap- 
pendix, the possibility of interpreting rigid spacetime 
symmetries as residual gauge transformations. Consider 
a free particle in in a gravitational background. The 
action is 



(Dl) 



What are the degrees of freedom? If we start with 
a generally covariant theory, we know that diffeomor- 
phisms correspond to gauge degrees of freedom. Let us 
consider the metric background as non-dynamical, take 
the passive view for the action of diffeomorphisms, and 
consider the simplified case where there exists a system 
of spacetime coordinates such that the metric is just 
Minkowski. Now we can make the following gauge fix- 
ing: We decide to stick with Minkowskian coordinates 
and only allow further diffeomorphisms if they keep this 
condition. We end up with Poincare transformations as 
residual gauge transformations. 

To count the degrees of freedom, we start with posi- 
tions and velocities i^, which amount to 8 degrees of 
freedom. To fix the r reparameterization invariance we 
must spend two constraints (for instance — t = and 
its stabilization i° = 1). We are left with 6 degrees of 
freedom. The residual gauge freedom consists of Lorentz 
transformations, but our gauge fixing forbids boosts and 
time translations, so we are left with three-translations 



and rotations. Rotations only affect 2 degrees of free- 
dom, because the norm of the velocity three vector is 
unchanged, and therefore we eliminate 5 of the 6 remain- 
ing degrees of freedom. We end up with a single degree 
of freedom that in our gauge fixing corresponds to the 
kinetic energy. 

Here we see a matter of interpretation as to what is 
and what is not a true degree of freedom. From our point 
of view, if the Lorentz invariance may be considered as 
the residual invariance found after a process of reducing 
the gauge group of general covariance, then the Lorentz 
degrees of freedom must be considered as gauge. But 
as residual gauge symmetries, they are not associated to 
constraints but to ordinary constants of motion. 

Summing up, the analysis of degrees of freedom, either 
in the case of Bianchi cosmologies or the simple case of 
a Minkowskian free particle, depends upon the point of 
view adopted. If one sticks to the formalism by itself, 
that is to say, to what the given Lagrangian yields as 
constraints, gauge transformations, and so on, one does 
not get the same number of degrees of freedom as if one 
considers that our theory comes from the reduction of 
an originally generally covariant theory. In this second 
point of view (which is the one we support) there are 
some symmetry transformations that are residual gauge 
transformations, that is, remnants of the reduction pro- 
cedure of the Lagrangian and the gauge group, which 
require gauge fixing. 
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